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Abstract: The ANTARES Collaboration is building a high-energy neutrino telescope at 2500 m depth 
in the Mediterranean Sea. The experiment aims to search for high-energy cosmic neutrinos through the 
detection of Cerenkov light induced by muons and showers resulting from neutrino interactions with the 
surrounding medium. The detector will consist of a three-dimensional array of 900 optical modules hous- 
ing photomultipliers. It will be composed of 12 strings, 5 of them being already in operation since January 
2007. The muon track is reconstructed from the arrival time and the charge of the signals obtained from 
the photomultipliers, whose positions are known by means of an acoustic positioning system. The recon- 
struction strategies include several steps among which there are: optical background filtering, algorithms 
for first estimations of the track parameters, and a final fit aiming to reach an angular resolution better 
than 0.3 degree above 10 TeV in the full detector Different reconstruction strategies will be presented 
and their application to the present real data analysis will be reviewed. 



Introduction 

When an ultra-relativistic particle (/3 ~ 1) moves 
in a medium, Cerenkov light is emitted at an angle 
depending on the refraction index n of the medium. 
In case of sea water n is about 1.34, and thus the 
Cerenkov emission angle becomes 

cos{ec) = - cosU2°). (1) 
n 

The emitted Cerenkov photons are detected by an 
array of photomultipliers installed at the bottom of 
the sea. Since January 2007 the ANTARES de- 
tector is a full three-dimensional array of photo- 
multipliers consisting of 5 strings detecting muons 
at a depth of 2475 meters below the sea level. 
ANTARES 10-inch photomultipliers are housed in 
pressure resistant glass spheres called optical mod- 
ules (OM). The 5 -line detector data taking has al- 
lowed the Collaboration to tune the various pro- 
cesses needed to collect data, and muon events 
have been detected on top of the optical back- 
ground (60-100 kHz most of the time). The cali- 
bration system has been proven to be efficient, and 
the knowledge of the charge, of the arrival time 
of the signals and the positioning of the OMs has 



allowed the first reconstruction codes to be tested 
in realistic conditions. The status of the experi- 
ment is discussed in [1]. Preliminary data studies 
show that a flux of atmospheric downward-going 
muons is triggering the detector at a rate of about 
1 Hz and that upgoing atmospheric neutrino can- 
didates have been identified. Among the various 
muon reconstruction algorithms under study, two 
different alternative methods are presented in this 
paper. Other well developed reconstruction meth- 
ods are described in [2], [3], and the implemen- 
tation of the simulated annealing algorithm [4] is 
under study. A discussion on event reconstruction 
techniques for Cerenkov neutrino telescopes can 
also be found in [5]. 

Track reconstruction strategies 

A minimum-ionising muon crossing the 
ANTARES detector causes the emission of 
Cerenkov photons, at an angle 9c of about 42° 
with respect to the track. As the muon energy 
increases, secondary particles are created by the 
muon energy loss along its path. These secon- 
daries produce photons which are delayed with 




Figure 1 : Schematic view of the kinematics related to 
the photon arrival time on the OM. The red line repre- 
sents the muon, while the blue line represents the pho- 
ton. 

respect to the light produced by the passage of the 
muon itself. Photons may also scatter in water, 
inducing an additional delay of the arrival times 
on the OMs. Assuming that the muon propagates 
as a straight line at speed c, the arrival time of the 
direct Cerenkov photon tf^, as shown in Fig. 1 can 
be directly linked to the muon track through the 
following equation: 

tth = t.,^Ui M I V M (2) 

* c \ tanOc J Vg \sin9c J 

where Vg is the group velocity of light, and £ and 
k are distances depending on the track parameters 
6, (j>, X, y, z, at t = tQ, where to is the time of the 
event. The track reconstruction methods usually 
consist in two main steps: a pre-fitting procedure 
where a first estimation of the track parameters is 
performed, and a final fit usually being a maxi- 
mum likelihood method, where a probability den- 
sity function (PDF) tries to include the full knowl- 
edge of the detector and the expected physics. This 
PDF is a function of the hit time residuals, de- 
fined as the arrival time differences between the 
observed hit times and the hit times expected for 
a direct Cerenkov photon travelhng directly from 
the muon to the OM 



U being the measured time, tj given by Eq. 2 for a 
given OM position Xi,yi, Zi, and a being the time 
resolution. A set of estimators, hit selections, fits 
and pre-fits is called a reconstruction strategy. 

Examples of reconstruction strategies 

Two strategies chosen as an illustration of the event 
reconstruction in ANTARES are presented here. 
The first strategy presented (here called strategy 1) 
is composed by a set of linear pre-fits, and each 
pre-fit is than taken as first guess for a final L1 — L2 
estimator fit as defined in [6]. The second strategy 
(here called strategy 2 or ScanFit) is composed by 
an iterative procedure giving a set of prefeiTed di- 
rections, each selected direction is then fitted with 
a minimum method. 

Pre-fitting procedures 

In general, the pre-fit step in track reconstruction 
must give an initial guess for the second minimi- 
sation step. The challenge is to find a relatively 
fast procedure giving a solution which is not too 
far away from the reality. Two pre-fitting proce- 
dures will be discussed in the following section: 
the linear and the scanning pre-fits. 

Linear pre-fit 

The line or Stenger fit algorithm [7] produces an 
initial estimation of the track parameters with the 
time and amplitude of the hits. This algorithm is 
efficient when the length of the track in the detector 
is large compared to the light attenuation length. 
The locations of each OM r^, recording a hit at 
time ti can be connected by a line 

r, ~ r + V (4) 

where r is the track position at < = and v the 
track velocity vector Under this assumption, a 
linear fit using space and time information of the 
event is performed through the minimisation of 

X' = ^a. (r. -r-vt,)' (5) 



(3) 



where Uh is the number of hits, is the measured 
charge for the hit i at position seen at time ti. 
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The factor Oj takes into account that a hit closer 
to the track will have a higher amplitude. The 
minimum is obtained from the parameter values r 
and V yielding vanishing partial derivatives. The 
resulting Eq. 5 can be solved analitically, giving: 



in) 



and 



where 



(r» • tj) - (rj) ■ {U) 



1 



nh 

\Xi) = / Qi Xi 

rih ^ 



(6) 



(7) 



(8) 



is the weighted mean of parameter Xi with respect 
to all hits. The Unear fit gives a vertex point r 
and a direction u = v/|v| at i = 0, where |v| 
is the mean speed of light propagating through 
the one-dimensional vector projection. Zenith 
and azimuth cf) angles of the fitted track are finally 
computed through the direction vector v and the 
(.T, y, z, 0,0, t = 0) parameters become the start- 
ing point for the subsequent minimisation step. 
The algorithm is applied to all hits in the event, and 
once the linear pre-fit solution is found, the track 
position is shifted to the plane perpendicular to the 
track and containing the barycentre of the detector 
giving a new position Fq . The time of the track is 
calculated accordingly: to = Tq /c. Then the av- 
erage t and the RMS of the hit time residuals with 
respect to the calculated track are evaluated. The 
average of this distribution is added to the track 
time in order to better match the Cerenkov model: 
^0 to +t- If the RMS of the time residuals distri- 
bution is more than 40 ns, the process is re-iterated 
and stopped when the condition is satisfied. The 
Unear fit gives a angular error of about 10° from 
the true solution for the 12- string detector. 

Scanning pre-fit 

Another possible way of finding the first guess for 
the track is an iterative search for the direction 6 
and (I) maximising the number of hits in the given 
direction which are close to the assumed track. 
Given the set of coincidences in the event, i.e. the 
hits i, j firing different OMs but the same floor and 
for which ti — tj is less than 20 ns, the set of hits 
composed by the earliest hits of the sample can be 



defined. For a given set of directions a ID clus- 
tering algorithm |8] i.e. a causality relation to de- 
termine space time correlations between the hits is 
applied to the set of earhest hits. Then, the same 
algorithm is applied a second time in order to add 
some other hits belonging to the event. The goal of 
this hit selection is to have a sample of hits having a 
high purity of direct Cerenkov light. The resulting 
set of hits is used to determine the three remain- 
ing coordinates x, y, z, using the full Hesse ma- 
trix. Only candidate directions with at least 5 hits 
are kept and passed on to the next stage, together 
with the selected hits. This procedure is more ac- 
curate than the Unear pre-fit and is capable of giv- 
ing a first guess which deviates only by about 5° 
from the true direction for the 12-string detector. 
The disadvantage is that this scanning procedure 
is computing time consuming. A similar iterative 
pre-fit was already discussed in [3]. 

Final fitting procedures 

From the first guess obtained by the pre-fit algo- 
rithm, a UkeUhood function £ is used to get a better 
estimate of the track parameters. This is normally 
done by the minimisation of —log{£). In general, 
the likeUhood function is the product of the PDFs 
Ci computed for each hit i 



(9) 



i=l 



The chaUenge is to find the best description of £j 
given an initial set of hits. Normally, the hit time 
residuals follow a complicated distribution, and 
a sophisticated parameterisation of the hit arrival 
time distribution is needed. Implementing a PDF 
in case of ANTARES requires a detailed knowl- 
edge of the detector response and of the photon 
propagation in water. The two strategies described 
don't make use of a specific PDF parameterisation 
but choose to use an approximated description of 
the PDF. The methods used to minimise are SIM- 
PLEX and MIGRAD for the first strategy and Ml- 
GRAD only for the second strategy (both can be 
found in the MINUIT package [9]). 



Fit with a LI L2 estimator 



The linear pre-fit can be relatively far from the so- 
lution. The efficiency of finding a good estimate 
with the final minimisation can be increased by 
rotating the direction and translating the position 
around the pre-fitted track. The first strategy pre- 
sented takes the set of linear pre-fits and the set of 
hits which fall one RMS from the initial track to 
begin the minimisation. The hit selection used is 
rather simple, so the problem is contaminated by 
the presence of a large number of outliers, i.e. a 
large number of scattered photons. In this particu- 
lar case, the solution can be found by minimising 
the LI — L2 estimator 

x2 = £2y'l + a.(tr)V2-2. (10) 
Fit with a estimator 

The second strategy presented takes the best solu- 
tions coming from the scanning pre-fit procedure 
and since the hit cleaning is particularly efficient, 
the hit time residuals already follow a Gauss distri- 
bution. In this particular case, the method of least 
squares (x^) can be applied to find of the best track 
parameters 

4 = 1 

As a result, the strategy gives a set of candidate 
fits. The candidate fits are ordered firstly on the 
number of hits and secondly on the normalised 
to the number of degrees of freedom. 

Comparison of the two methods 

The two different strategies were applied to the 
same set of ANTARES 5-line data. Each event has 
been carefully calibrated in time, charge and posi- 
tion. The first strategy presented has a reconstruc- 
tion efficiency on triggered events of 85% while 
the second one of 91%. The zenith angle distribu- 
tions resulting from the two reconstruction meth- 
ods are shown in Fig. 2. In green the distribu- 
tion obtained for the first strategy (linear pre-fit 
and L 1 -L2 estimator fit) and in red the distribution 
obtained for the second strategy (scanning pre-fit 








Figure 2: Zenith angle distributions given in arbitrary 
units resulting from the two reconstruction strategies de- 
scribed in the paper. In green the distribution obtained 
for the first strategy presented (linear pre-fit and L1-L2 
estimator fit) and in red the distribution obtained for the 
second strategy presented (scanning pre-fit and fit)- 
No analysis cuts are applied. 

and x^ fit). No analysis cuts were applied. The 
two different methods presented do not make use 
of any detailed likelihood description, they use ap- 
proximations with analytical estimators, but these 
preliminary results obtained are encouraging. The 
two strategies can be improved with the inclusion 
of a more detailed description of the PDF. 
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